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We calculate the photon production cross section arising from the hard scattering of partons in 
nucleon-nucleon collisions by taking into account the intrinsic parton transverse momentum distri- 
bution and the next-to-leading-order contributions. As first pointed out by Owens, the inclusion 
of the intrinsic transverse momentum distribution of partons leads to an enhancement of photon 
production cross section in the region of photon transverse momenta of a few GeV/c for nucleon- 
nucleon collisions at a center-of-mass energy of a few tens of GeV. The enhancement increases as 
■y/s decreases. Such an enhancement is an important consideration in the region of photon momenta 
under investigation in high-energy heavy-ion collisions. 
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I. INTRODUCTION 

Currently, high-energy heavy-ion collisions are used to 
produce matter under extreme conditions. One of the ob- 
jectives is to search for the quark-gluon plasma (QGP) 
which is expected to exist at high temperatures and/or 
high baryon densities. (For an introduction to this field, 
see @-§|.) Photons arising from the electromagnetic in- 
teractions of the constituents of the plasma will provide 
information on the properties of the plasma at the time 
of their production. Since photons are hardly absorbed 
by the produced medium, they form a relatively 'clean' 
probe to study the state of the quark-gluon plasma. The 
detection of these photons during the quark-gluon plasma 
phase will be of great interest in probing the state of the 
quark-gluon plasma, if it is ever produced |fj|-|l5|] . 

Photons are also produced by many other processes in 
heavy- ion reactions. They can originate from the decay 
of 7T° and if. As n° particles are copiously produced in 
strong interactions between nucleons, photons originat- 
ing from the decay of w° are much more abundant than 
photons produced by electromagnetic interactions of the 
constituents of the quark-gluon plasma. The photons 
from the decay of tt° and 77° can be subtracted out by 
making a direct measurement of their yield, obtained by 
combining pairs of photons. Because of the large num- 
ber of 7T° produced, this subtraction is a laborious task, 
but much progress has been made to provide meaningful 
results after the subtraction of the photons from the ir° 
and rp backgrounds Jll] [f|. Photon measurements ob- 
tained after the subtraction of the photons from meson 
decays are conveniently called measurements of "direct 
photons" . 

Direct photons are produced from the interaction of 



matter in the QGP phase, a mixed QGP and hadron 
phase, a pure hadron gas, and hard QCD processes PB. 
Different processes give rise to photons in different mo- 
mentum regions. We are interested in photons from the 
quark-gluon plasma which are found predominantly in 
the region of photons with low transverse momentum, 
extending to the region of intermediate transverse mo- 
mentum of 1-3 GeV/c. At photon energies up to 2 GeV, 
photons can come from the decay of 7r° and rf resonances 
as well as from p, u>, rj\ and a\ , and the interaction of the 
hadron matter via np — > 7/5 and 7T7t — > reactions (see 
for example [p|-|T^]). One may wish to go to the region 
of photon transverse momentum p-yx > 2 GeV/c. If a 
hot quark-gluon plasma is formed initially, clear signals 
of photons from the plasma could be visible by examining 
photons with j? 7 t in the range 2 — 3 GeV/c ||^-[l0|]. On 
the other hand, photons in this region of transverse mo- 
menta are also produced by the collision of partons of the 
projectile nucleons with partons of the target nucleons. 
Such a contribution must be subtracted in order to infer 
the net photons from the quark-gluon plasma sources. 

Recent investigations on photon production by parton 
collisions in hadron-hadron collisions include the work by 
Aurenche, Baier, Douiri, Fontannaz, and Schiff [^8|Jl9[ 
and by Baer, Ohnemus, and Owens |20(| , who have per- 
formed QCD calculations up to second order in a. Ex- 
tensive comparisons with experimental data have been 
carried out covering a large range of incident energies 
and photon transverse momenta. In these QCD calcula- 
tions, as well as other similar QCD calculations the 
intrinsic transverse momenta of the partons have been 
neglected. 

Because of the finite size of the transverse dimension 
of a nucleon, one expects that the partons in a nucleon 
have an intrinsic transverse momentum distribution with 
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a width of the order of y/2h/0.5 fm. The intrinsic trans- 
verse momentum distribution affects the distribution of 
the produced photons. Previous investigations of the 
photon transverse momentum distribution in leading- 
order (LO) calculations with a constant i-T-factor already 
indicate the importance of the intrinsic transverse mo- 
mentum of the partons p^ ]. For p 7 x from 3 to 8 GeV/c, 
the photon production cross sections for hadron-hadron 
collisions at at 19.4 GeV < y^s < 63 GeV are enhanced 
when the intrinsic transverse momentum distribution is 
taken into account, in agreement with earlier measure- 
ments. These findings are further supported by recent 
experimental and theoretical investigations for hadron- 
proton collisions at yfs =31.7 and 38.8 GeV H). The 
enhancement increases as y/s decreases. Therefore, in 
the region of our interest involving photons with p y T of 
about 2-3 GeV/c at a center-of-mass energy of about 20 
GeV, as in the measurements of the WA80 and the WA98 
Collaborations j[l],[l2), the intrinsic transverse momen- 
tum of partons plays an important role and cannot be 
neglected. How the photon transverse momentum dis- 
tributions at different y/s are affected by the parton in- 
trinsic transverse momentum is the main subject of our 
investigation. 

There is also another important effect which is impor- 
tant in photon production. Previous calculations of pho- 
ton production indicate the importance of higher-order 
QCD terms jl8|-^0|]. The next-to-leading-order calcula- 
tions lead to a correction factor, the if -factor, with a 
magnitude of about 2. It is necessary to include the ef- 
fects of the next-to- leading-order corrections. One ex- 
pects that the next-to-leading-order effects are nearly in- 
dependent of the intrinsic transverse momentum. One 
can include the effects of the intrinsic transverse momen- 
tum and the next-to-leading-order corrections by first 
treating them separately and then combining them as 
independent multiplicative factors. The final results can 
be compared with experimental data. Good agreement 
with experimental data at different reaction energies will 
form the basis of extrapolation to the region of interest 
in high-energy heavy-ion collisions. 

Previous work by Cleymans et al. [pl| on photon pro- 
duction in pp collisions covers a region from y/s = 23 
GeV to 1.8 TeV. These authors use the next-to- leading- 
order calculation program of Aurenche et al fl8|| , without 
assuming parton intrinsic transverse momentum. How- 
ever, the results of these authors cannot be interpreted 
as the conclusive evidence for the absence of any parton 
intrinsic momentum effect at energies of y/s ~ 20 GeV, 
because their theoretical results are lower than the UA6 
data at y/s — 24.3 GeV by about a factor of 2 (see Figure 
3 of UW), and the theoretical results are lower than the 
E706 data at y/s = 30.6 GeV by about a factor of 1.7 
(see Figure 4 of pl|). As mentioned above, Owens ]2^ ] 
has earlier explained similar discrepancies in the photon 
production data of pp, n + p, and ir~p reactions at 19.4 



GeV< y/s < 63 GeV as arising from the parton intrin- 
sic momentum p2| . We shall provide further support of 
Owens' observations by analytical and numerical studies. 

In this manuscript, we first examine photon produc- 
tion in nucleon-nucleon collisions and extend our consid- 
erations to nucleus-nucleus collisions. In nucleon-nucleus 
and nucleus- nucleus collisions, there are shadowing ef- 
fects for photon production in nuclei. Previous exper- 
imental investigations of photon production using ir~ 
beams on nuclei suggest that the effective shadowing in 
the production of photons can be represented by a target 
mass dependence of the form A a with a ~ 1.0 p5[ , which 
can be explained theoretically j26) . Because a is close to 
unity, the effective shadowing for photon production in 
nuclei is quite weak. We shall not include the shadowing 
effect on photon production in nucleus- nucleus collisions. 



II. PHOTON PRODUCTION IN THE 
HARD-SCATTERING MODEL 

In the relativistic hard-scattering process, the produc- 
tion of a photon in a hadron-hadron reaction (h a + ht — > 
7 + X) arises from the direct interaction of a parton of 
the hadron h a with a parton of the other hadron h^. The 
basic processes are the Compton process g + q (or q) — > 
7 + q (or q), and the annihilation processes qq — ► jg 
til- 

One represents the probability for finding a parton a 
with momentum a=(x a , clt) in the hadron h a by the par- 
ton distribution function G a /h a (x a , ax), which depends 
also on the momentum transfer Q 2 . The x a and Q 2 de- 
pendences of the parton distribution function Gw/, a for 
various partons have been obtained by fitting large sets of 
experimental data in different hard-scattering processes 

In leading-order perturbative QCD, the cross section 
for h a + hi, — > 7 + X is a convolution of the parton mo- 
mentum distributions and their elementary collision cross 
sections E 7 d 3 a(a b — > jX')/dp 3 



E. 



/ dx a da,Tdxi,dbT 



d 3 a(h a h b -> jX) 
dpi 

' LLU 

n i \n i u \c d 3 a{ab -> jX') 
xG a / fta (x Q ,a T )G b/hi; (a; b ,b T )£; 7 — . (1) 



The parton invariant cross section is related to da/dt by 



E 



d 3 a(ab -> 7V) da(ab -> jX') 



dp. 3 



7 



ivdt 



5(s + i+u-^2m 2 ) 



i=l 



s_V ~ ( m i + ™2) 2 }{s - (mi - m 2 ) 2 } , 
[s + ra\ — m\) 

where mi, m,2 are the masses of the incident partons, 
and m.3, and 7714 are the masses of the outgoing particles. 
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In the above equation, we have used the Mandelstam 
variables 



s = (a + 6) 2 , 
t= (a-p 7 ) 2 , 
u = (b-p^) 2 . 



(3) 



For dcfinitcness, we shall use the nucleon-nuclcon 
center-of-mass system to refer to the momenta of the 
photon and the partons. After averaging the flavors, the 
colors, and the isospins of the initial states of q, q, and 
g, and summing over final states, the differential cross 
section da(a + b — > 7 + X')/dt are 

daigq — ► 7$ _ 1 fe q \ 2 8na s a e 



dt 



6 \e J (s — to 2 ) 2 
,2 x 1 



s — m 

„2 



u — m 
™2 



4 V u — m 2 



,(4) 



and 



d<r(qq — > jg) _ 4/e g \ 2 87ra 5 a e f 

di ~ 9\e) s{s- 4m 2 ) ]\t - m 2 



t — to 



m 



u — TO' 
„2 



4 \ u — m 2 



t — m 2 



(5) 



where to is the quark mass, e q and e are the quark charge 
and the proton charge respectively, and a e and a s are the 
electromagnetic fine-structure constant and the strong 
interaction coupling constant, respectively. The strong 
coupling constant a s is related to the momentum transfer 
Q 2 by 



a s (Q 2 ) 



12tt 



(33-2Ar / )ln(Q 2 /A 2 )' 



(6) 



where Nf is the number of flavors, and A is the QCD 
scale parameter. 

We shall assume a factorizable parton distribution 
function where the intrinsic transverse momentum dis- 
tribution can be factored out in the form 



G a /h a (x a ,a T ) = F a/ha (x a ,Q 2 )D a (a T ). 



(7) 



For convenience, we normalize the transverse momentum 
distribution D a (a,T) such that 



J da T D a {a T ) = 1. 



(8) 



Then, the function F a /} la (x a , Q 2 ) is the usual parton dis- 
tribution function without assuming a parton intrinsic 
transverse momentum distribution, as given by HHo). 



III. EFFECT OF PARTON INTRINSIC 
TRANSVERSE MOMENTUM 

Quarks, antiquarks, and gluons are asymptotically free 
at high momenta, but their properties at low momenta 



are governed by the confinement of these particles inside 
a hadron. As a consequence, the transverse momentum 
distribution of these particles is given by their momen- 
tum wave functions. The width of the momentum distri- 
bution is related to the inverse of the size of the radius 
of confinement. Thus, if one takes a confinement radius 
of 0.5 fm, the standard deviation a of the parton intrin- 
sic momentum should be of the order of Ti/ (0.5 fm)=0.4 
GeV/c, and the width of the intrinsic transverse momen- 
tum about a/2 x 0.4GeV/c=0.65 GeV/c. 

In the usual PQCD calculations, one neglects the par- 
ton intrinsic transverse momentum by approximating the 
parton transverse momentum distribution D(ax) with 
a delta function D(ay) = 5(ay). Then, the integration 
over ax and 6<r in Eq. (Q) can be trivially carried out. 
In such an approximate description, the colliding partons 
have zero transverse momentum. The produced photon 
acquires a transverse momentum through the basic colli- 
sion processes of gq(oT q) — ► 7?(or q) and qq — * jg, with 
differential distributions as given by Eqs. (||) and (||). 

The transverse momentum of the photons produced 
from the collisions of the constituents of the quark-gluon 
plasma expected in high-energy heavy-ion collisions lies 
predominantly in the low p 7 T region, extending to the re- 
gion of intermediate p^T of a few GeV/c. Hard-scattering 
processes in nucleon-nucleon collisions also produce pho- 
tons in this intermediate transverse momentum region. 
It is necessary to determine the hard-scattering contri- 
butions to the photon production in order to extract in- 
formation on the quark-gluon plasma photons. The de- 
termination of the hard-scattering photons will require 
the inclusion of the effects of the intrinsic parton trans- 
verse momentum. 

Previously, the effect of parton intrinsic transverse mo- 
mentum has been investigated by Owens [^2| in hadron- 
nucleon collisions at different energies. In the range of 
3 GeV/c < p 7 T < 8 GeV/c, Owens found by numeri- 
cal calculations that the parton intrinsic px leads to an 
enhancement of the photon production cross section and 
the enhancement increases as the incident hadron energy 
decreases [P2fl . The enhancement of photon production 
cross section due to the parton intrinsic transverse mo- 
mentum has also been reported in recent numerical cal- 
culations for hadron-proton collisions at -y/s =31.7 and 
38.8 GeV [p4| . It is instructive to show analytically in 
this section how the intrinsic transverse momentum en- 
hances the photon production cross section. We shall 
make simplifying assumptions in this section in order to 
carry out the six-dimensional integration in Eq. (^) ana- 
lytically. 

We consider the production of a photon with y 7 = 
and transverse momentum p~ ( x- We introduce the photon 
light-cone variable, x 7 = (E 1 + p 7Z )/y/s. At y 1 = 0, we 
have Py Z = and x 7 = p-yT / \fs- The photon momentum 
can be represented by 
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7 = { E -y, P^t, V~iz) = {x-yy/s, p jT , 0) . 



(9) The corresponding x a is given from Eq. (Raj by 



The momenta of the target parton a and projectile parton 
b can be written as (see Eq. (4.10) and (4.11) of |§) 



a = \ x 



2 1x a \fs 



, it, - x a - 



y/s a, 



2 2x a y/s 



y/s b% y/s 

Xb— + ^ j=, b T , x b — - -= 

2 2xb\/s 



b\ 



2 2xb\fs 
The Mandelstam variable s, t, and u are then 



, CLrpbrp 

S = X a XbS H 2clt • by-, 



x b H ) x^s + 2b T ■ P 7 t, 

x s 



( 4\ 
t = — I x a H I ir 7 s + 2a^ • P 7 t • 



(11) 

(12) 
(13) 
(14) 



We take the partons to be massless. The momentum 
conservation condition, s + t + u — 0, leads to a relation 
of x a in terms of x b given by 

a% 1 



Xb + {b^/xbs) — x 1 



2 2(aT + br) ■ p 7 T — 2aT ■ br x a b\ Xba\ 



+ - 



XbS X a S 



To simplify the last two terms in the square bracket, we 
note that x a ~ (see Eq. ( P8|) below) . The above equa- 
tion becomes 



|p 7 T - a T 



6t| 2 



{x& + (b^/xbs) — x 7 }4 



(15) 



which can be solved for x a to obtain x a — x a {xb). We 
assume a parton distribution function of the form F(x) = 
A(l — x) n . The integration over the light-cone variables 
x a and Xb in Eq. (Fy) becomes 



/ 



dx a dx b A a A b (l - x a ) n "(l - x b ) 



„ b S (jg J(Xq - XgjXb)) 

tt dt s{xb — X~/) 
da 



TTs(Xb — Xj) dt 

where 

f(x b ) = n a \n{l -x a {x b )} +n b ln(l - x b ). 
The function f(xb) is an extremum of Xb at 



x b 



XbS 



l n a (l - x b ) |p 7 t - Qt - &r| 
n fc (l - a; a ) v/s 



(16) 



(17) 



(18) 



X a + 



I n b (l - x a ) \p~ (T - Qt - &r| 

7l a (l-Xi,) y/s 



(19) 



(10) The integration over x a and in Eq. (|T^) becomes 



dx a dx fc ^ a A(l - x a ) na {l - x b ) nb -^-5(§ + i+u) 

tt at 

l2ir A a A b S \n h d(J /on\ 

-{1-XaT) na {l-X bT ) ""-J-, (20) 



B ITS {x b T — Xj) 



dt 



where x b r and x a r are the solution of x b and x a from 
solving Eqs. (18) and ( |l9| ) respectively, and B is 



_9V 
(9x^ 

b x b =x bT 

n b 2^/n a n b (l - 2x 7 ) + x 7 (n a + n b ) 
(1 - 2x 1 ) 2 



(21) 



Consider the case with p 7 x >> c, where a is the standard 
deviation of the parton px distribution. We can expand 
in powers of the intrinsic transverse momenta to get 



\p lT -a T -b T \ ~p 7 t(1 



p jT ■ (gT + b T ) 



(22) 



Neglecting terms of second power in ay and br and using 
Eqs. (fs|), (|l^), and (p2|), we can separate out x a T and 
x b T into 



£iT = x i0 + A 4 



(t = a, b), 



where 



XaO 



l n b (l - x a0 ) 
n a (l - Xbo) 



f 71 b (l - x a0 ) p lT ■ {a-T + b T ) 
n a (l-x bQ ) p~fT\fs 



n a (l - x b0 ) 

Xbo=[l + \ —r. r X 7 , 

V n b (l - Xao) 



and 



' n a (l - x b0 ) p 7 t ■ jar + b T ) 
n b (l - x a0 ) p^ry/s 



(23) 



(24) 



(25) 



(26) 



(27) 



For the Compton process (gq — > jq), n a ~ 6 and 7i& ~ 4 
p7{ , and for the annihilation process (qq — > 73), n a = 
n;, ~ 4. Therefore we have approximately 



Xao ~ a^fco ~ 2x 7 . 



(28) 
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We are now ready to carry out the integration over the 
intrinsic transverse momenta ax and for in Eq. ([!]). The 
dependence of (1 — x a x) na on the intrinsic momentum 
or can be separated out to be 



(l-x aT ) n ° ~ (1 - x a0 ) n ° exp{- 



n a A a 
l — x„ 



-} 



~ (1 - a*,)- expfv^ y (a T + 6 ^ -} (29) 
I (1 - x a0 )p lTX /s J 

The factor (1— x b T) nb can be similarly separated. The in- 
tegral over the transverse momentum in Eq. (Ill) becomes 



Jda T db T D a (a T )D b (b T ){l - x aT ) n *(l - x bT ) r 
= (1 - 2x 7 )"«+"" J da T db T D a {a T )D b {b T ) 



i <i<7 



xeM2V^ ?7 T \ (aT , +bT l ] d ^ (30) 



(1 — 2x~ ( )p~ fTy /s J (it 

where we have used Eq. (^8|) . We assume a Gausssian 
parton transverse momomcntum distribution, D(j>t) = 
exp{— py/2cr 2 }/27rcr 2 with <j?y>= 2<r 2 , and we obtain 



fda T db T D a (a T )D b (b T ) expi2y/n a n b f 1 ^ ("^ + M 
J L {yfs-2p lT )p lT ) d< 

y da T db T D a (a T -\)D b (b T - A)^, 

(31) 



= exp 



4cr 2 n a n fc 
(v^-2p 7 t) 2 



where 



_ 2a 2 y /n a n b p lT 

(y/s - 2p lT )p jT ' 

The first factor on the right-hand side of Eq. ( |3l|) , 



f 4cr 2 n Q n b 
cxp<^ - 



(32) 



(33) 



is the enhancement factor due to the parton intrinsic 
transverse momentum and the variation of the parton 
distribution function. We can understand the origin of 
this enhancement in the following way. The integration 
over x a and x b leads to a cross section proportional to 
(1 - x aT ) n "(l - x bT ) nb (Eq. ©). If there were no par- 
ton intrinsic transverse momentum, x a T would be given 
by x a o ■ The presence of the intrinsic transverse momen- 
tum leads to x a x varying about x a o with the variation 
A a depending on the intrinsic transverse momentum (Eq. 
(p3|)). Because of the steep decrease of (1— x a ) Ua (l—x b ) rib 
as a function of x a and x b , the averaging over the trans- 
verse momenta leads to the enhancement factor K\. 

The enhancement factor K\ arising from the parton 
distribution function allows us to determine the scale of 
p-yT at which the effect of the parton intrinsic pr will 
be important. The enhancement factor k\ is close to 



unity if 2 < p\ > n a n b « (y/s — 2jj 7 t) 2 and is much 
greater than unity if 2 < p\ > n a n b is large or compa- 
rable to (y/s — 2p lT ) 2 . As n a n b ~ 24 for the Comp- 
ton process, the importance of the intrinsic px depends 
on whether 1 \J <p\ > is small compared to \fs — 2p~ / x- 
When 1\J <p\ > is large or comparable to (y/s — 2p 1 r), 
the enhancement factor k\ becomes quite large. This 
is the main origin for the large enhancement factor for 
p 7 T < 6 GeV/c in hadron-proton collisions at y/s ~ 20 
GeV. 

For the cross section da/dt, the dominant contribu- 
tions come from terms proportional to 1/t and 1/u. It 
can be approximated by 



£~c(i + I 

dt \t u 



From Eqs. @, (|T|), and @, we have 



and 



i ~ |p 7 t ~ b T \ 2 , 



Wit - a T \ 



(34) 



(35) 



(36) 



The integration over the intrinsic transverse momentum 
distributions can be carried out to yield 



[da T db T D a (a T — \)D b (b T — \)C ( — ; 

J \(p 7T -a T 



1 



) 2 (p 



It 



P 



It 



2C 



n 2 
P 7 T 



{p lT -X) 2 -2a 2 
where first factor on the right hand side 



K'2 



(p 7 T-A) 2 -2 ( r 2 ' 



(37) 
(38) 



is the enhancement factor due to the parton intrinsic 
transverse momentum and the variation of da/dt. This 
enhancement factor «2 arising from da / dt is close to unity 
\ip lT »a. From Eqs. (p3J), (|35|), and (|36|), the enhance- 
ment factor due to da/dt becomes much greater than 
unity when p^ is of the order \J <p\ >. 
Putting all the results together, we have 



E. 



d 3 a(h a h b -> >yX) 



2tt AA a A bPlT 




B TT^i 
where n is given by 

Aa 2 n a n b 



1 



n a +n b 



2C 

2 " 

PjT 



k = K1K2 = exp 



(y/7s-2 PlT ) 2 j p 2 T (l-X) 



r> 2 

P-yT 



2a 2 



(39) 



(40) 



The quantity k in Eqs. ( p9| ) and ([4C]) is the ratio of the 
photon production cross section with the parton intrinsic 
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transverse momentum to the the cross section without 
the parton intrinsic transverse momentum. From Eq. 
( ^p| ) we infer that forj^T >> \J < p\ > , which allows 
the expansion in Eq. n is greater than unity and 

represents an enhancement due to the parton intrinsic 
transverse momentum. It is a function of both -^/i and 
p 7 T- It decreases as y/s increases. For the same y/s, the 
first factor increases with increasing p 7 T while the second 
factor increases with decreasing p 7 r ■ 

If one uses 2cr 2 = 0.9 (GeV/c) 2 as given by Owens 
then for p lT = 6 GeV/c at yfs = 20 GeV, Eqs. Q and 
( |38| ) give K\ = 2.2 and n 2 — 1.27, and a combined en- 
hancement factor k — 2.8, in rough agreement with the 
enhancement factor of about 3.8 obtained from numerical 
calculations in Figure 3. There is a substantial enhance- 
ment of the photon production cross section at yfs ~ 20 
GeV due to the parton intrinsic transverse momentum. 
For p 7 T — 6 GeV/c at yfs = 63 GeV, n decreases to 1.07, 
in agreement with the numerical results in Figure 4. This 
value of k is only slightly greater than unity, indicating 
that the parton intrinsic pt has only small influence for 
p lT ~ 6 GeV/c at sfs > 63 GeV. 

The enhancement factor k in Eq. (|4^) shows that one 
should expect a large enhancement due to the intrinsic 
transverse momentum for the production of photons with 
transverse momentum p 7 r close to y/s/2 such that y/s — 
2p*yT is not much greater than y / <p|3>. In this case, 
x a o r~j 2p~ / T/yJ r s is close to unity and the magnitude of 
the cross section is however quite small. 

The above analytical integration of the six-dimensional 
integral of Eq. ([j]) serves the purpose of illustrating the 
important features of the enhancement due to the par- 
ton intrinsic transverse momentum. In order to compare 
with experimental data quantitatively, we need to carry 
out numerical integrations of Eq. ([!]) in their full com- 
plexities without simplifying assumptions in the follow- 
ing sections. 



IV. TRANSVERSE MOMENTUM 
DISTRIBUTION OF PARTONS 

The momentum distribution of partons for small px 
depends on confinement, while the momentum distribu- 
tion at high px depends on the constituent counting rule 
which is a power law in nature |u 32 ||. They have dif- 
ferent functional forms. To carry out our investigation 



of the parton intrinsic momentum in numerical calcula- 
tions, we join the two different functional forms at ptm 
and assume the parton px distribution function of the 
form 



D(pr) = C 



1 



1 + e( pT ~ p ™)/ A 
1 

1 g(pTJV/-PT)/A 



D^Pt) 
D 2 (pr) 



(41) 



where C ~ 1 will be determined by the normalization 
condition Eq. (||). 

Following Owens [^2| , we use the Gaussian distribution 
D\(pt) for the distribution at low px as given by 



-Di(Pt) 



1 



7T(p T 2 ) 



-Pt/<Pt> 



(42) 



For large pr, we take the transverse momentum distri- 
bution D 2 (px) to be given by a power-law distribution 
as in the spectator counting rule of Blankcnbecler and 
Brodsky |3| 



D 2 (pt) 



(Pt 2 



M 2 y 



(43) 



where a = D 1 {p T m){Ptm + M 2 ) n . The results of our 
calculation for photon production are rather insensitive 
to this transverse distribution D 2 . Using information 
from the work of Sivers, Blankenbecler, and Brodsky |32| , 
we choose n = 4 and M — 1 GeV to describe the trans- 
verse momentum distribution D 2 - For numerical pur- 
poses, we set the matching transverse momentum to be 
Ptm = 2 y/ < p T > , and the matching width A = 0.25 
GeV/c. 



V. NEXT-TO-LEADING-ORDER 
PERTURB ATIVE QCD 

It is well known that higher-order perturbative QCD 
corrections are important in the evaluation of the cross 
sections for many hard-scattering processes (see for ex- 
ample, j3^j2^]). For photon production, contributions 
up to the second order in a s have been considered by 
Aurenche, Baier, Douiri, Fontannaz, and Schiff, and by 
Baer, Ohncmus, and Owens ||l8|-po||. These calculations 
involve the evaluation of contributions from all next- 
to-leading-order Feynman diagrams. In these calcula- 
tions [fl8|-po[ , the intrinsic transverse momentum of the 
partons has been neglected by using a delta function 
distribution <5(pt) so that integrations over the trans- 
verse momentum of the partons can be trivially carried 
out. From these leading-order plus next-to- leading-order 
(LO + NLO) calculations, the effects of the higher-order 
terms can be represented by the if-factor defined by 



hii>.) = \E^(LO- 



NLO) 



(44) 



where E 1 (d 3 a/dp'^)(LO + NLO) is the photon obtained 
with the inclusion of both the lowest-order and the next- 
to-leading-order terms, and E 1 (d 3 a/dp^)(LO) is the pho- 
ton cross section obtained in the lowest-order calculation. 

Numerically, the if-factor is about 2 to 3.5 and is 
a slowly varying function of the photon transverse mo- 
mentum (see Figs. 1 and 2). It represents an effective 
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modification of the strength of the coupling constant at 
the vertex, as in a vertex correction, due to the initial- 
state interaction between partons [j33|-[36"|] . On the other 
hand, the vertex correction depends on the relative mo- 
mentum between the colliding partons |34j-E3(|, which in 
turn depends mainly on the longitudinal momenta; the 
vertex correction is insensitive to the parton transverse 
momenta. Because of this property, it is reasonable to 
consider the A-factor to be approximately independent 
of the parton transverse momentum; that is 



d 3 a 

E 1 — (LO + NLO - 



rf 3 



PT) 



' E 1 ^r {LO + NLO) 



E 1 — {LO+NLO) I E^ — {LO) 



1 dp% 



(45) 



where the abbreviation "PT" stands for calculations 
where the parton transverse momentum distribution is 
taken into account. Upon using the above approxima- 
tion to include the effects of the next-to-leading-order 
corrections and the parton intrinsic transverse momen- 
tum distribution, the photon production cross section is 
given by 

' ^ NLO + PT)^K{p 1 )E 1 —r{LO+PT). (46) 

dp 6 
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Fig. 1. The if -factor as a function of the photon 
transverse momentum p 7 r at y 1 — 0. (a) for pp 
collisions at y/s = 17.3 GeV and (b) for ^/s = 19.4 
GeV. 



VI. PHOTON PRODUCTION IN 
NUCLEON-NUCLEON COLLISIONS 

From Eq. (f46|), we can take into account the effects of 
intrinsic transverse momentum and next-to-leading-order 



contributions by carrying out the calculation in two steps. 
We perform the LO+NLO calculations, without the in- 
trinsic parton transverse momentum distribution, using 
the numerical code from Aurenche et al. 18 ljj. The 
ratio of the LO + NLO calculation to the leading-order 
(LO) calculation then gives the A-factor as a function of 
the photon momentum, as defined by Eq. ((44|). We then 
carry out separately a numerical calculation of the lead- 
ing order with the intrinsic transverse momentum dis- 
tribution. Equation ( pH ) is then used to give the cross 
section when both effects are included. 

The photon cross section depends on the functional 
dependence of Q 2 in terms of physical quantities such as 
p 7 T of the photon. Owens has studied different choices of 
the functional form, from Q 2 = p 2 T to Q 2 = p 2 T /4, and 
found that Q 2 — p 2 T /2 gives the best description of the 
experimental data for production of photons with large 
transverse momenta [^2|. We shall accordingly use this 
functional dependence for our investigation. 

We first carry out the LO and the LO+NLO calcula- 
tions using three different parton distribution functions: 
(1) the Duke and Owens parton distribution function (set 
1) \27\jj2) the CTEQ4 parton distribution function (set 
3) 2jj]29||, and (3) the MRS96 parton distribution func- 
tion (set 1) J30|. The ratio of the LO+NLO result to the 
LO result gives the if -factors which are shown in Figs. 
1 and 2, and the the results for LO and LO+NLO cross 
sections are shown in subsequent figures. The A"-factor 
in Figs. 1 and 2 as a function of the photon transverse 
momentum p 7 T for y 7 = indicate that it has a value 
of approximately 2 for 1 GeV/c < p 7 T < 6 GeV/c in 
nucleon-nucleon collisions at y/s — 17.3 and 19.4 GeV. 
It behaves slightly differently for y/s = 63 GeV, where it 
drops from about 3 to 1.4, as p y T decreases from 2 GeV/c 
to 6 GeV/c. It becomes nearly a constant for p^r > 6 
GeV/c. 




DO 

1.5 f- CTEQ4 

MRS96 

1.0 



2.0 



4.0 



6.0 8.0 

Pir (GeV/c) 



10.0 



12.0 



Fig. 2. The A-factor as a function of the photon 
transverse momentum p y T at y 7 = for pp collisions 
at V« = 63 GeV. 

We next carry out a lowest-order calculation with the 
inclusion of the parton transverse momentum distribu- 
tion of the partons. Previously, Owens noted that the 
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transverse momentum distribution of dilepton pairs gives 
a direct measure of the parton intrinsic transverse mo- 
mentum [p2[ . Data obtained by Kaplan et al. |37j for pp 
collisions at ^/s — 27. A GeV indicate < p^{n + H~) >= 
1.9 (GeV/c) 2 which corresponds to a parton transverse 
momentum < p\ >= 0.95 (GeV/c) 2 in the lowest or- 
der approximation. The value of < p\ > decreases if 
higher-order terms are included p3fl . The convolution 
method of Q is however a model-dependent prescrip- 
tion for incorporating the effects of the intrinsic trans- 
verse momentum. Nevertheless, it indicates that the 
amount of intrinsic transverse momentum needed to de- 
scribe the data is influenced by the amount of QCD dy- 
namics included in the analysis. The analysis of Owens 
using the lowest-order QCD with a constant iT-factor 
indicates that <p\>= 0.9 (GeV/c) 2 gives a good agree- 
ment with photon production data in pp, ir + p and ir^p 
collisions at 19.4 GeV< y/s <62 GeV {22j. Recent anal- 
ysis Q suggests even a greater value of <p\ > of about 
2 (GeV/c) 2 . More future work is needed to clarify the 
situation on the the magnitude of the parton intrinsic 
transverse momentum. 



t 

f s" 2 = 19.4 GeV, y=0 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 



PyT (GeV/c) 

Fig. 3. The photon invariant cross section for 
nucleon-nucleon collisions at ^fs — 19.4 GeV and 
y-y — calculated with different <p\ > values of par- 
ton intrinsic momentum distribution and the Duke 
and Owens parton distribution function (set 1). The 
NLO corrections have been included. 

To see the dependence of the photon transverse mo- 
mentum distribution on the magnitude of < pip > of the 
parton intrinsic transverse momentum, we show in Fig. 



3 the results for pp collisions at y/s — 19.4 GeV, using 
the Duke and Owens parton distribution function (set 
1). The cross sections for <p\ >= are less than those 
for < p\ >= 0.9 (GeV/c) 2 by about a factor of 3 to 4, 
which are less than those for <p\ >=1.8 (GeV/c) 2 by 
about a factor of 2 (except at p 7 T — 1 GeV/c). There 
is a substantial enhancement of the cross section due 
to the parton intrinsic transverse momenta. The value 
of < p\ >= 0.9 (GeV/c) 2 provides a good description 
of the experimental data. The enhancement factor for 
<Pj, >= 0.9 (GeV/c) 2 is relatively constant for y/s ~ 20 
GeV. 




0.0 2.0 4.0 6.0 8.0 10.0 12.0 



p jT (GeV/c) 

Fig. 4. The photon invariant cross section for 
nucleon-nucleon collisions at y/s = 63 GeV and 
y 7 = calculated with different <p\ > values of par- 
ton intrinsic momentum distribution and the Duke 
and Owens parton distribution function (set 1). The 
NLO corrections have been included. 

The enhancement of the photon production cross sec- 
tion in the region of p^T under consideration poses an in- 
teresting puzzle. If the running coupling constant could 
be taken as a fixed constant, the total cross section ob- 
tained by integrating the photon momentum should be 
independent of the width (or <p\ >) of the parton intrin- 
sic transverse momentum distribution. An enhancement 
of the cross sections in all regions of p 7 T appears to vi- 
olate this normalization of the total photon production 
cross section. The puzzle is resolved by noting that the 
cross section is not enhanced in all regions of p 7 t- The 
photon production cross section is enhanced in the re- 
gion p 7 T > > v <Pt >, and it is suppressed for p^r < < 



8 



\J <p\ > because partons acquire energies through their 
intrinsic px and the probability for producing very low 
energy photons therefore decreases when partons have 
intrinsic px- A hint of the suppression for low j» 7 t can 
be found in Figs. 3 and 4 where for p^r < 1-2 GeV/c the 
cross section for the case <p T >— 1.8 (GeV/c) 2 is lower 
than that for <p\ >= 0.9 (GeV/c) 2 . While we expect 
a suppression for this region of small p lT « <p\ >, 
it is a different matter whether perturbative QCD can 
remain valid for this region where the associated value of 
Q 2 becomes very small. We shall not consider this region 
of small p-yT- 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 



p v (GeV/c) 

Fig. 5. The photon invariant cross section for 
nucleon-nucleon collisions at ^fs — 19.4 GeV and 

2/7 = °- 

Figure 4 shows the direct photon production results for 
nucleon-nucleon collisions at ^/s — 63 GeV and y 1 = 0. 
The nucleon-nucleon collision data from R806/807 @ 
and R108 (CCOR) Collaborations ||| arc shown in the 
figure. For p 7 T > 1-2 GeV/c, the cross section is en- 
hanced as < p T > increases. The magnitude of the en- 
hancement is smaller than that for y/s = 19.4 GeV, indi- 
cating that the effect of intrinsic transverse momentum 
decreases with increasing y/s. For p 7 x < 1-2 GeV/c, the 
cross section initially increases as < p\ > increases but 
decreases as <p T > increases further, indicating that the 
trend of enhancement is reversed for p 7 x « \J <Pr >■ 
The value of <p T >= 0.9 (GeV/c) 2 provides a good de- 
scription of the experimental data. We shall henceforth 
follow Owens to use <p^>— 0.9 (GeV/c) 2 when we 



include the parton intrinsic transverse momentum distri- 
bution in our subsequent calculations. 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 



Pir (GeV/c) 

Fig. 6. The photon invariant cross section for 
nucleon-nucleon collisions at y/s — 17.3 GeV and 
y 1 = 0. The scale on the right refers to the photon 
distribution per central Pb-Pb collision event in the 
WA98 measurement Q . 

In Fig. 5 we present results of the photon cross section 
for the LO calculations, the LO+NLO calculations, and 
the LO+NLO+PT results including further the effect of 
the transverse momentum distribution, for pp collisions 
at ^ = 19.4 GeV and y 1 = 0. The LO + NLO calcula- 
tions are enhanced from the LO calculations by about a 
factor of 2, and the intrinsic momentum effect brings an- 
other factor of about 4 to 8, depending on the parton dis- 
tribution function. The intrinsic transverse momentum 
of partons enhances the cross sections by a substantial 
factor and is very important in the region of intermediate 
Pjt. The data for ,/s = 19.4 GeV from E704 Q, E629 
|9|, and NA3 ^ are also shown in Figs. 3 and 5. (The 
NA3 data are taken from the reaction p + 12 C — > 7 + X 
where we assume a nuclear dependence of A , as in 
Ref. @.) The WA80 data in Figs. 3 and 5 represent 
the upper limits of the cross section. The experimental 
data in the region of p 7 T from 2 to 4 GeV/c can be de- 
scribed well when both the effects of higher-order terms 
and the parton intrinsic transverse momentum are taken 
into account. For the region p^x > 4 GeV/c, the experi- 
mental data of NA3 and E629 are different. The data of 
NA3 above 4 GeV/c fall within the predictions using the 
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CTEQ4 and the MRS96 parton distribution functions. 
Within the considerable experimental uncertainties, the 
experimental data agree with the LO+NLO+PT predic- 
tions, although there are small differences in the results 
from different parton distribution functions. 

Recent photon production experiments using high- 
energy heavy ions were performed for Pb-Pb collisions 
at 158 GeV per nucleon |I^,[l3|], which correspond to 
nucleon-nucleon collisions at y/s — 17.3 GeV. In order 
to facilitate the comparison with experimental data, we 
show the theoretical results for pp collisions at yfs = 17.3 
GeV and y 7 = in Fig. 6. The intrinsic transverse mo- 
mentum enhances the cross sections by a factor of about 
4 to 8. It is therefore an important effect which must be 
included in the discussion of photon cross sections in the 
region of p^r between 1 and 6 GeV/c. 



VII. PHOTON PRODUCTION FROM THE QCD 
HARD-SCATTERING PROCESSES IN 
NUCLEUS-NUCLEUS COLLISIONS 

We can determine the direct photon production cross 
section in high-energy nucleus-nucleus collisions arising 
from the hard scattering between partons. Because of 
the weak effective shadowing effect in photon production 
by parton collisions as indicated by experimental pA data 
p5| , we shall neglect the shadowing effect. The photon 
distribution in a collision of a nucleus A and a nucleus B 
at an impact parameter b due to the hard scattering of 
partons is given by 



E,^^ -(b) = n NN (b) E 



dN^ N {H.S.) 
d Pl y ~> '" NNy "' ~ 7 d^ 



(47) 



where n NN (b) is the number of inelastic nucleon-nucleon 
collisions, 



n NN (b)=AB T AB {b)a. 



NN 
in 1 



(48) 



and the photon distribution in an inelastic nucleon- 
nucleon collision is 



dN^ N {H.S. 

dp-y 



— Ery 



da" N {H.S.) 



r N N 



dp. 



(49) 



In the above equations, a^ n N and a^ N are the nucleon- 
nucleon total inelastic cross section and photon produc- 
tion cross section respectively, T AB (b) is the thickness 
function for the nucleus-nucleus collision 101 



T AB (b) = / db A db B T A {b A )T B {b B )t(b -b A b B ), 



(50) 

where T A (b) and T B (b) are respectively the thickness 
functions for nucleus A and B, and the baryon-baryon 



thickness function t(b — b A — b B ) can be approximated 
by a delta function, t(b — b A — b B ) = 5(b — b A — b B ). 
Note that when one combines Eqs. (|47|), (|8|), and (|I|), 
the nucleon-nucleon inelastic cross section a^ n N cancels 
out and one obtains the result 



dN AB 



dp 



(b) = ABT AB (b)E 7 — 2— _1 L. (51) 



dp. 



Writing the result for nucleus- nucleus collisions in the 
form of Eq. ( f47| ) provides a more intuitive picture of the 
scaling of the photon multiplicity from nucleon-nucleon 
collisions to nucleus- nucleus collisions. 

Experimentally, one measures the yield of photons per 
event of nucleus-nucleus collisions in coincidence with the 
measurement of the degree of inelasticity of the events, 
as characterized by the multiplicity of produced parti- 
cles. Results are given in terms of the photon yield per 
inelastic collision event for the most-inelastic fraction / 
of all inelastic events. The most-inelastic fraction / de- 
pends on the maximum impact parameter b m , and one 
can deduce a relation of / as a function of the maximum 
impact parameter b m from the Glauber model by 



f(bm) 



Cdb{l-(l-T AB (b)ag N ) AB } 
f o °°db{l-(l-T AB (b)a£ N )^} 



(52) 



For measurements within the most-inelastic fraction 
f(b m ) of events, the photon distribution per inelastic 
event is given by taking the average of Eq. ([l7|) over the 
range of impact parameters from to b m . We there- 
fore have the photon distribution per inelastic event in a 
nucleus-nucleus collision given by 



Ey 



dN* B (H.S) 
dp 1 



{b)=<n NN > (b m )E. 



da" N {H.S) 



7 (j NN 



dp-. 



(53) 



where the average number of nucleon-nucleon collisions 
per event within the range of impact parameters from 
to b m is 



<n >(b)- 1^ dbAB T AB (b)ag- 
NN tdb { l- { l-T AB[ b)^VB } 



(54) 



In Fig. 7(a) we show < njvw > (b m ), the number 
of nucleon-nucleon inelastic collisions averaged over the 
range < b < b m , as a function of the maximum impact 
parameter b m for Pb-Pb collisions, and in Fig. 7(b) the 
corresponding most-inelastic fraction f(b m ). 

The photon production data from the WA98 Collab- 
oration were collected for the most-inelastic 10% of the 
total number of Pb-Pb inelastic collisions |l^] . For this 
value of f(b m ) = 0.1, the impact parameter b m is found 
from Fig. 7 to be 4.44 fm. With an inelastic nucleon- 
nucleon cross section of 29.4 nib 0, the average number 
of inelastic nucleon-nucleon collisions within this range 
of impact parameters from to 4.44 fm is <hnn >=654. 
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This is the scaling number which one must multiply to 
the photon yield in a nucleon-nucleon inelastic collision 
to give the total photon yield per central Pb-Pb collision. 




K (fin) 

Fig. 7. (a) The number of inelastic nucleon-nucleon 
collisions in Pb-Pb collisions averaged over the im- 
pact parameter from to b rn as a function of the 
impact parameter b rn . (b) The fraction of the total 
inelastic cross section for collisions with an impact 
parameter from to b m . 

Using this scaling factor, we obtain the photon distri- 
bution from hard scattering of partons per central Pb-Pb 
collision at a nucleon-nucleon collision energy of y/s — 
17.3 GeV and y 1 = 0. They are shown as the dashed and 
dashed-dot curve in Fig. 8 for the CTEQ4 parton distri- 
bution function (set 3) and the MRS96 parton distribu- 
tion function (set 1). They have been calculated for the 
most-inelastic 10% of the inelastic Pb-Pb collisions. They 
should be compared with photon distributions from the 
quark-gluon plasma at different temperatures to be dis- 
cussed in the next section and shown as the solid curves 
in Fig. 8. 

The scaling between nucleon-nucleon collisions to 
nucleus-nucleus collisions can be extended for S-Au col- 
lisions which were studied by the WA80 Collaboration 
0. The WA80 photon data were collected for the 7.4% 
of the most-inelastic events for which the maximum im- 
pact parameter b m is found to be 2.9 fm and the number 
of nucleon-nucleon inelastic collisions <tinn> averaged 
over the range of impact parameters from to b m is 160. 
This scaling number has been used to convert the WA80 



upper limit E^dN^ Au /d 3 p 7 data to E^da™ /d 3 p 7 data 
in Figs. 3 and 5. 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 

p v (GeV/c) 



Fig. 8. The invariant photon distribution per cen- 
tral Pb-Pb collision event in the WA98 measurement 
p2[ at y 1 — as a function of p 7 T- The dashed and 
dashed-dot curves are results from parton hard scat- 
tering obtained in LO+NLO+PT calculations, and 
the solid curves are from a quark-gluon plasma with 
different initial temperatures To and a critical tem- 
perature T c = 200 MeV. 



VIII. PHOTON PRODUCTION FROM THE QGP 
IN NUCLEUS-NUCLEUS COLLISIONS 

In a quark-gluon plasma, the constituents of the 
plasma interact with each other. The interaction between 
the constituents leads to the production of photons via 
gq(ov q) — > 7q(or q) and qq — > 75 reactions, which are 
the same processes as in photon production in parton 
collisions |6]- |l0| , |l4| ] . Upon taking the lowest-order Feyn- 
man diagrams in the evaluation of the cross section and 
assuming a thermal mass mth = gT/V6, the rate of pho- 
ton production in a quark-gluon plasma at temperature 
T is § [Eq. (16.66) of §] 

dA 7 5a e a s 2 J3.7388.EU 

= 9^ fMT ln ri^-f ' (55) 

where g 2 /Air — a s is the strong- interaction coupling con- 
stant, and fq(py) is the Fermi-Dirac distribution of the 
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quarks in the plasma which, for E 7 >> T, is approxi- 
mately the Boltzman distribution 



fq(Pl 



-E^/T 



One can introduce the A'-factor to take into account 
the next-to-leading-order effect. The X(_E 7 )-factor as 
shown in Fig. 1(a) obtained for parton collisions (at 
^/s =17.3 GeV) can be used as an approximate if-factor 
for constituent collisions in the quark-gluon plasma. The 
photon distribution per event is then 



E~, 



dN 1 _ 5 K(E 1 )a e a s 
~ ~ 9 



dp 



2ir 2 



d A xf q {p- 



)T 2 In 



3.7388K 



9 2 T 



(56) 



We shall assume Bjorken hydrodynamics in which the 
proper time varies with temperature T as 



T 



(57) 



where tq is the initial proper time when the system is at 
a temperature To. We consider a system with an initial 
volume of V given by Atq where A is the average over- 
lap area for the range of impact parameters considered. 
Then, during the evolution of the system from temper- 
ature To to the critical phase-transition temperature T c 
in the quark-gluon plasma phase, the number of photons 
emitted can be obtained by carrying out the integration 
in Eq. @. The result is 

E dNy = 5g(^)a e a s ( VtqT^ 1 



dp y 67r 2 \ h 4 

£i(z) + e- z ln 



To E y 



3.7388z 
Anas 



z—E^/Tq 



z=E^/T a 



(58) 



where E\(z) is the exponential integral J z °° dte~ t /t fl43| . 

The photon distribution from the quark-gluon plasma 
for T c =200 MeV and different values of To is shown in 
Fig. 8 for the 10% most-inelastic central collisions of Pb 
on Pb, where we have taken the initial time parame- 
ter to to be 1 fm/c. As one observes, the photon num- 
bers increase as initial temperature increases. This arises 
because it takes a longer period of time to cool a sys- 
tem to the phase transition temperature if the system 
is initially at a higher temperature. The photon num- 
ber distributions from hard scattering of partons are also 
shown in Fig. 8. One observes that the photon distri- 
bution from the quark-gluon plasma is greater than that 
from the hard scattering when the initial temperature of 
the plasma exceeds about 0.28 GeV. 
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Fig. 9. The total invariant photon distribution per 
central Pb-Pb collision event in the WA98 measure- 
ment at y 7 = as a function of p 7 t, including the 
contributions from nuclcon-nucleon hard scatterings 
and the quark-gluon plasma with an initial temper- 
ature To. 

In Fig. 9, we show the prediction for photon produc- 
tion when both the hard-scattering contributions and the 
quark-gluon plasma contributions are added together, for 
the 10% most-inelastic central collisions of Pb on Pb at 
158 GeV per nucleon. The hard-scattering results are ob- 
tained from the LO + NLO + PT calculations. The total 
distribution depends on the plasma initial temperature. 
The have been calculated by assuming a critical temper- 
ature T c = 200 MeV. They may be used to compare with 
experimental data. 



IX. DISCUSSIONS AND CONCLUSIONS 

In order to detect photons which are emitted during 
the quark-gluon plasma phase, it is necessary to deter- 
mine the contribution of photons by non-quark-gluon 
plasma sources. The hard scattering between partons 
also produces photons in the momentum region of inter- 
est. It is therefore important to study the hard-scattering 
contribution to photon production in the intermediate 
p 7 T region of a few GeV/c. 

For photon production at *Js « 63 GeV and p 7 r 
up to about 6 GeV/c, the effects of the intrinsic trans- 
verse momenta of the partons cannot be neglected. We 
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have examined how the intrinsic transverse momentum 
distribution of partons affects the transverse momentum 
distribution of the photons produced in nucleon-nucleon 
collisions. The intrinsic momentum of the partons leads 
to an enhancement of the cross section by a factor rang- 
ing from 3 to 8 for photons with intermediate p-yT- Such 
an effect is an important consideration in photon mea- 
surements under investigation in high-energy heavy-ion 
collisions. 

We have also compared the magnitude of photons pro- 
duced by hard scattering with photons produced by the 
quark-gluon plasma. We find that for photon production 
in the region of 1-3 GeV/c in Pb-Pb central collisions at 
158 GeV per nucleon the quark-gluon plasma contribu- 
tion exceeds the hard-scattering contribution when the 
initial temperature is greater than about 0.28 GeV. 
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